Heavy flavor nuclear modification factor: more baryons than mesons less energy loss 
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The suppression of the nuclear modification factor for heavy flavor hadrons is usually attributed 
to the energy loss of heavy quarks propagating in a QCD plasma. Nevertheless it is puzzling that the 
suppression is as strong as for light flavors. We show that when accounting for the quark momentum 
shift associated to the opening of the recombination/coalescence channel for hadron production in 
the plasma, it is not necessary to invoke such strong energy loss. This shift is expressed in terms of 
an increase of the heavy baryon to meson ratio in nuclear with respect to proton collisions. When 
this mechanism is included along with a moderate energy loss, data from RHIC and LHC for the 
nuclear modification factor of electrons coming from heavy flavor decays as well as of charm mesons, 
can be reasonably described. 



PACS numbers: 25.75.-q, 25.75.Nq 

The suppression of single hadron transverse spectra in 
nuclear collisions, with respect to a superposition of inde- 
pendent proton collisions, is one of the main results from 
the BNL Relativistic Heavy Ion Collider (RHIC) and the 
CERN Large Hadron Collider (LHC). This suppression 
is quantified in terms of the nuclear modification factor 
(Raa) and one of its main features is that heavy flavor 
hadrons are equally suppressed as light hadrons [J, 
Such behavior was first obtained from the analysis of elec- 
trons from the decay of heavy flavors and later confirmed 
from the analysis of charm mesons [H, Q . 

When the suppression of heavy flavors is only at- 
tributed to energy loss in the QCD medium, the above 
result is surprising for if the main contribution comes 
from radiative processes, the dead cone effect [j| should 
prevent heavy quarks from losing as much energy as light 
ones. This motivated the reviewing of energy loss sce- 
narios to incorporate contributions from collisional pro- 
cesses, diffusion, geometry, as well as dynamical proper- 
ties of the medium Q. However even these refined sce- 
narios do not yet provide a fully convincing explanation 
for the properties of the heavy flavor Raa- 

Much less attention has been paid to the fact that a 
shift of the hadron momentum in the nuclear medium can 
come not only from a loss of energy but also from a mo- 
mentum redistribution when the quarks from the medium 
form either baryons or mesons. This is the central idea 
behind the recombination/coalescence scenario as a new 
channel for hadron production in a heavy-ion environ- 
ment Q • In average, the three quarks forming the baryon 
come from lower momentum bins than the two quarks 
making up a meson. Since there are more quarks with 
lower momenta there is a larger chance to form baryons 
than mesons. Transverse flow increases the effect since 
this makes the momentum distribution for heavier par- 
ticles (hadrons) to fall less steeply than for lighter ones 



(mesons). A direct consequence of this momentum redis- 
tribution is an increase of the baryon to meson ratio in 
nuclear with respect to proton collisions. This ratio has 
been measured for a large variety of light and strange 
hadrons in high-energy nuclear collisions @ . The upshot 
is that for intermediate transverse momenta, the ratio is 
enhanced with respect to the corresponding one in pro- 
ton collisions. Although no measurements exist for the 
case of heavy flavors, there are model calculations that 
describe this enhancement 0-[H| . To test this scenario in 
quantitative terms, we use one of these models, the Dy- 
namical Quark Recombination Model (DQRM) [12], to 
compute the heavy baryon to meson ratio which in turn 
is used to compute R e AA an d Raa- We show that when 
this increase is accounted for, only a moderate energy 
loss is needed to reproduce the data. 

For definitiveness, let us concentrate on describing the 
nuclear modification factor for a single heavy flavor, say 
charm (c) quarks. The number of c-quarks produced in 
nuclear (AA) or proton (pp) collisions in a given momen- 
tum bin can be obtained from counting the correspond- 
ing number of hadrons with c-quarks, namely, the num- 
ber of open charm mesons (N AA , ), charm baryons 

( N aa I PP ) and h idden charm mesons (N A C A / pp ) 
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Normalizing the proton case to the average number of 
binary collisions (n&) we have that, accounting also for 
a possible shift in energy e, the number of produced c- 
quarks in one and the other environments, must satisfy 
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that is 



«4 + N^ A + Nf A ) = e(n b ) (N? p + + Ng). (3) 
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^From Eq. ([3]), wc can build the nuclear modification 
factor for _D-mesons in terms of the number of charm 
mesons and baryons in pp and AA collisions and get 
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The last term in Eq. (j4]) can be written as 
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Using Eq. ([5]) into Eq. (j4]), we can write 
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the 



where we have defined rj 
";i tin of hidden charm to 

Npp/N^, is very small, to an excellent approximation 
we can rewrite Eq. ([6]) as 
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Therefore, even in the absence of energy loss (e = 1) the 
nuclear modification factor for D mesons is smaller than 
one, provided the ratio of charm baryons to open charm 
mesons is enhanced in AA with respect to pp collisions. 

The same enhancement is responsible for the suppres- 
sion of the nuclear modification factor for heavy-flavor 
electrons. For definitiveness, let us again focus on elec- 
trons originating from the decay of charm quarks. In a 
given momentum bin R e AA can be expressed as [ToL fl3| 
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where B D ' A ^ e is the branching ratio for the decay of D 
mesons and charm baryons into electrons, respectively. 
Using Eq.([7|), we can write Eq. ([5]) as 
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FIG. 1: (Color on line) DQRM charm baryon to meson ra- 
tio in AA compared to the same ratio in pp collisions. The 
DQRM curves are computed for two transverse expansion ve- 
locities v t — and 0.65. For the ratio in pp we use PYTHIA 
simulations at ^/s]7n = 200 GeV and 2.7 TeV with 35 x 10 6 
and 15 x 10 6 events, respectively. Shown are also fits to the 
simulations. 



where, in order to introduce the energy loss factor e, wc 
have used Eq. ([3]) ignoring the contribution from hidden 
charm mesons. Also x = B A ^ e / ' B D ^ e and the function 



(7) T% A is given by 
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It has been shown that when x < 1, T AA is also 
smaller than one when the ratio of charm hadrons to 
open charm mesons is enhanced in AA with respect to 
pp collisions [l(|. Therefore Eq. © states that even in 
the absence of energy loss, the nuclear modification factor 
for single electrons is smaller than one, provided the ratio 
of charm hadrons to open charm mesons is enhanced in 
AA with respect to pp collisions and that electrons are 
more copiously produced from open charm mesons than 
baryons (a; < 1), which is indeed the case. 

In the DQRM the probability V to recombine quarks 
into mesons and baryons depends on density and tem- 
perature and thus on the proper time r describing the 
evolution of the heavy-ion reaction up to hadronization. 
The evolving probability differs for hadrons made up by 
two and three constituent quarks and is computed by a 
variational Monte Carlo simulation. The relative popula- 
tion of one or the other kind of cluster at low densities can 
be fixed by combinatorial arguments (see Refs. [l(| [H| 
for details). This model is well suited to describe baryon 
and meson production and its ratio at low and intermedi- 
ate pt- The hadron transverse momentum distribution in 
central AA, assuming Bjorkcn dynamics and transverse 



3 



DQRM e =0.55 
DQRM E =0.4 
D" 
D* 




4 5 6 
p(GeV) 



10 



DQRM 2.7 TeV e =0.55 
DQRM 2.7 TeV e =0.4 
DQRM 0.2 TeV £ =0.55 
DQRM 0.2 TeV e =0.4 
STAR 
ALICE 




FIG. 2: (Color on line) Nuclear modification factor for D 
mesons compared to ALICE data. The curves are computed 
using the DQRM with vt = 0.65 as appropriate for LHC en- 
ergies. The charm baryon to meson ratio in pp is taken from 
the PYTHIA simulation described in Fig. [1] For simplicity, 
the energy loss parameter is taken as two constant values, 
e — 0.55 (upper curve) and e = 0.4 (lower curve). 

velocity expansion v t , is given by 

dN m t Ay Pnncl p 

— ; — r~ = 9 — a a — / t(xt P[t\ 

Ptdptdy 47T At J tq 

x / (p tS inh7 7t /T) e - m ' cosh '"/ T , (11) 

where m t is the transverse mass, Ay the rapidity interval, 
|0 nucl the nuclear radius, At = tj — To the proper time 
interval and T the proper time dependent temperature 

2 

T = T o(^)', (12) 

with Vg = 1/3. Jo is a Bessel function I of order zero. 
Vt and rjt are related through vt = tanhr^. g is the 
degeneracy factor that takes care of the spin degree of 
freedom. 

Fig. [JJ shows the DQRM charm baryon to meson ra- 
tio. We set the masses of the charm baryon and mesons 
to m A = 2.29 GeV and m D = 1.87 GcV. We take the 
initial hadronization proper time tq = 1 fm, at an initial 
temperature To = 175 MeV and the final hadronization 
temperature 7/ = 100 MeV that, according to Eq. (JT2"|) . 
corresponds to 17 = 8 fm. Shown are the cases between 
Vt = and v t = 0.65. The figure shows also the baryon 
to meson ratio in pp at ^/snn = 200 GcV and 2.7 TeV, 
obtained from PYTHIA simulations with 35 x 10 6 and 
15 x 10 6 events, respectively. Shown are also fits to the 
simulations. Notice that, as expected, the charm baryon 
to meson ratio is enhanced in AA with respect to pp col- 
lisions. 

Fig. [U shows R% A compared to ALICE data [|. The 
theoretical curves are computed using Eq. (UJ) with the 
heavy baryon to meson ratio obtained in ^L4 from the 
DQRM with the same parameters as before and the par- 
ticular value vt — 0.65, which is a standard choice for 



FIG. 3: (Color on line) Nuclear modification factor for non- 
photonic single electrons compared to STAR and ALICE 
data. The curves are computed using the DQRM with 
vt = 0.65 for the LHC case and v t = 0.55 for the RHIC 
case. The charm baryon to meson ratio in pp is taken from 
PYTHIA simulations for ^/s]7n = 200 GeV and ^/s^T =2.7 
TeV for the sum of charm and beauty hadrons. For simplic- 
ity the energy loss parameter is taken as two constant values, 
£ = 0.55 (upper curves) and e = 0.4 (lower curves). 

the transverse expansion velocity at LHC energies. The 
heavy baryon to meson ratio in pp is obtained from the 
PYTHIA simulation shown in Fig. [1] for LHC energies. 
To see the effect of the energy loss parameter, for simplic- 
ity, we take two constant values, s = 0.55 (upper curve) 
and e = 0.4 (lower curve). We notice that even in this 
simple scenario, data are well described and the energy 
loss parameter does not need to be as small as in the case 
of light flavors, which in this language means e ~ 2, to 
account for the suppression in R AA . 

Fig. [3] shows R AA compared to data from STAR 
and ALICE Q. The theoretical curves are computed 
using Eq. (0) with the heavy baryon to meson ratio 
obtained in AA from the DQRM and in pp from the 
PYTHIA simulations of Fig. [TJ To account for the find- 
ing that electrons from heavy flavor decays come al- 
most in equal proportions from the decays of charm and 
beauty hadrons for p t > 5 GeV [3], here we consider 
a single species of heavy baryons and mesons with ef- 
fective masses. We take mP = 3.57 GeV, the aver- 
age between the masses of the D° and the B° mesons 
and m A = 3.95 GeV, the average between the masses 
of the A c and the A^. Also, we consider that the 
possible charm and beauty mesons decaying into elec- 
trons or positrons are D ± (B D± ^ e± = 16%), D°, D° 
(B D °> D °^ e± = 6.53%), Df {B D ^ e± = 8%) and B ± 
(B B± ^ e± = 10.8%), B°, B° {B B ° • B °^ e± = 10.1%). 
The possible charm and beauty baryons decaying into 
electrons or positrons are A c , A c (B Ac ' A <= _! ' e = 4.5%), 
At and Af, (B Ab ' A >>^ e = 5.35%) (the experimentaly re- 
ported branching ratio coreresponds to the scmilcptonic 
decay A^ —> A c I £7. Here we consider that half of this 
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comes from the decay into electrons). The other param- 
eters used for the AA case are as before, with v t = 0.55 
for the RHIC case and Vt = 0.65 for the LHC case. For 
simplicity, the energy loss parameter is taken also as two 
constant values e = 0.55 (upper curves) and e = 0.4 
(lower curves). Once again, even in this simple scenario, 
data are well described for p t > 2 GeV and the energy 
loss parameter does not need to be as small as in the case 
of light flavors to account for the supression in R AA . For 
p t < 2 GeV, the rise in the data is usually attributed 
to other effects like shadowing, which is not considered 
in our approach. The model curves are not significantly 
affected if the effective masses are slightly varied. 

In conclusion, we have shown that when accounting 
for the medium's quark momentum redistribution when 
these recombine/coalesce to form mesons and baryons, 
the heavy flavor nuclear modification factors can be de- 
scribed without the need of a large energy loss. This 
momentum redistribution is encoded in the increase of 
the heavy baryon to meson ratio. We emphasize that the 



results are valid, independent of the model as long as the 
baryon to meson ratio increases in AA with respect to 
pp collisons. This increase is expected based on general 
grounds, since it represents a feature of the openning of 
the recombination/coalescence hadron formation channel 
in AA collisions. Upcoming upgrades to RHIC detectors 
and to ALICE are expected to increase the capability to 
directly look at this quantity and thus experimental tests 
of the mechanism advocated in this work will be available 
in the near future. 
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